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We have performed ab-initio calculations of the electronic structure and exchange couplings in the 
layered vanadates CaV^Os and MgV2 0s. Based on our results we provide a possible explanation 
of the unusual magnetic properties of these materials, in particular the large difference in the spin 
gap between CaV2 0s and MgV2 0s. 



I. INTRODUCTION 

Spin- 1/2 ladder models can describe the magnetic be- 
havior of a variety of quasi-one-dimensional systems [ 1 1 . 
Examples include the cuprate materials SrCu203 [2|, 
LaCuC>2.5 H and (Sr,Ca)i4Cu2404i Q. Spin excitations 
in the isolated ladders have a finite energy gap, which 
makes them prototype spin liquids. This is of interest 
in relation to high temperature superconductivity, since 
upon doping they become resonating-valence-bond liq- 
uids, with a spin excitation gap and a dominant quasi- 
long range pairing correlations Q . 

Another example of the spin- 1/2 ladder systems are 
the layered vanadate compounds CaV^C^ and MgV^Os. 
Although CaV^Os and MgV^Os have nearly identical 
vanadium oxygen planes however their magnetic prop- 
erties are strikingly different. CaV^Os has a large spin 
gap of about 600 K ||, while the spin gap in MgV^Os 
is very small only about 20 K ||. In contrast to the 
planar cuprates, where a hole in the Cu x 2 — ?/ 2 -orbitals 
results in a strong antiferromagnetic exchange coupling 
for the 180° bonds and a weak ferromagnetic one for the 
90° bonds, the exchange interactions in these vanadates 
can be more complicated as shown in Fig. 1. Even the 
signs of the many exchange couplings are not obvious in 
these materials. So one has to resort to ab-initio numer- 
ical calculations to get information about the relative as 
well as absolute values of the exchange couplings in these 
systems. The determination of the exchange couplings is 
crucial to understand the markedly different spin gap be- 
havior in these compounds. 

In this paper, we shall report on the ab-initio calcula- 
tion of the exchange couplings using the LDA+U method 
and discover that they are indeed different in these two 
compounds, consistent with their magnetic properties. 
As the various exchange couplings are related to the bare 
hopping matrix elements, we shall extract them using 
a recently developed systematic downfolding scheme . 
The advantage of the downfolding method is that only 
the important orbitals referred to as the active chan- 



nels are retained in the basis and the rest are down- 
folded thereby providing a single or few band tight bind- 
ing model capable of reproducing the details of the LDA 
bands close to a prescribed energy, which is usually the 
Fermi energy. We shall use such few band models to 
extract the various hopping matrix elements. This will 
form a basis to understand the widely different spin gap 
behavior in CaV^Os and MgV^Os. 

The remainder of the paper is organized as follows: In 
section II we shall briefly recapitulate the crystal struc- 
ture of CaV^Os and MgV^Os. In section III we shall 
present our LDA+U calculations for the various exchange 
couplings and compare them with available experimental 
results. Section IV will be devoted to electronic struc- 
ture calculations based on the TB-LMTO method [||, 
followed by the downfolding method to extract the var- 
ious hopping matrix elements, in order to explain the 
different exchange couplings in these materials. Finally 
the conclusions are given in section V. 




FIG. 1. Various magnetic interactions in Ca(Mg)V20s 
compounds. Two nearest vanadium ladders with different 
z-coordinates are shown in different shades. Jl, the ex- 
change interaction between nearest V atoms is ferromagnetic 
for CaV205 and antiferromagnetic for IVhgV^Os. J2 and J3 
are respectively the antiferromagnetic exchange interactions 
along the rung and leg of the ladder. J4 is the antiferromag- 
netic exchange interaction between the V atoms along the 
diagonal of the ladder. The magnetic structure used in the 
calculations is marked with up- and down- arrows. Half of 
Jl and all of J4 exchange interactions are frustrated in this 
magnetic structure. 
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FIG. 2. On the left the crystal structure of CaV2 0s and 
on the right the crystal structure of MgV^Os. Ca and Mg 
atoms are not shown, large balls represent V and small - O 
atoms. In each figure the oxygen atoms constitute a pyramid 
which are linked by edge and corner sharing as described in 
the text. 

II. STRUCTURE 

The main building block of the crystal structures of 
Ca(Mg)V2 05 compounds are the V ions roughly in the 
center of a pyramid of oxygen ions as can be seen in 
Fig. 2. The crystal structure of CaV20s is primitive 
orthorhombic with space group Pmmn and lattice con- 
stants a=11.35 A, 6=3.60 A, and c=4.89 A. As shown in 
Fig. 2 (left), the structure is formed by a linkage of VO5 
pyramids having apex oxygens in the direction of the c- 
axis. Oxygen edge- and corner-shared zigzag V chains 
are formed along the 6-axis, where the nearest neighbor 
V-V distance is 3.03 A. Along the a-axis, these chains 
are linked by sharing corners with the V-V distance of 
3.49 A. It thus forms a quasi two dimensional ladder layer 
in the ab plane with the leg along the zigzag V chains 
(i.e along 6) while the rung in the perpendicular direc- 
tion (i.e. along a). The Ca atoms are located between 
the layers and are surrounded by eight O atoms. 

The crystal structure of MgV205 is base centered or- 
thorhombic with space group Cracm; and lattice con- 
stants a=11.02 A, 6=3.69 A, and c=9.97 A. Again, the 
structure can be described as a linkage of VO5 pyramids 
having apex oxygens in the direction of the c-axis as can 
be seen in Fig. 2 (right). The V zigzag chains extend along 
the a-axis by sharing edges and corners of the pyramids 
and the nearest neighbor V-V distance is 2.98 A. They 




FIG. 3. A pair of corner shared VO5 pyramids. The various 
angles indicated in the figure are quoted in Table 1. 

are also linked by sharing corners along the 6-axis, with 
the V-V distance of 3.37 A, which again leads to a quasi- 
two dimensional ladder layers in the ab plane. However in 
contrast to CaV205, these layers stack alternately with 
the distance of (c/2) and the Mg atoms are located be- 
tween the layers and each are surrounded by six oxygen 
atoms. As a consequence, there is a puckering of the 
V2O5 layers to accommodate Mg ions in the tetrahcdral 
co-ordination. 

All the structural data that we have discussed above 
are summarized in Table 1 and indicated in Fig. 3. It 
should be noted that as a consequence of the puckering 
of the V2O5 layers the tilting angle of the corner shared 
pyramids a (see Fig. 3) is appreciably smaller in MgV205 
in comparison to CaV20s. 



III. THE LDA+U AND EXCHANGE 
COUPLINGS. 

The LDA+U method was shown to give good re- 
sults for insulating transition metal oxides with a par- 
tially filled d-shell || . The exchange interaction parame- 
ters can be calculated using a procedure based on the 
Green function method, developed by Lichtenstein et. 
al. ]TT| , p~2[ . This method has been successfully applied 
to calculate the exchange couplings in KCUF3 jl2) and in 
layered cuprates jl3| . 

The LDA+U method §,|l§ is the Local Density Ap- 
proximation (LDA) modified by a potential correction 
restoring a proper description of the Coulomb interac- 
tion between the localized d-electrons of transition metal 
ions. This is written in the form of a projection operator: 



TABLE I. Structural data for CaV20s and MgV20s compounds. The lattice constants (in A) are a along the rung of the 
ladder, 6 - along the leg, and c - in the vertical direction. The distances (in A) between nearest V ions between ladders, along 
the leg, the rung and the diagonal are denoted as d„„, d; e9 , A rung and ddi ag , respectively. The angle a (Figure 3) is between 
V-O-V where oxygen atom is placed between two vanadium atoms forming the rung. The angle f3 is that between O-V-O 
where one is rung oxygen and the other is apical oxygen. Finally the angle 7 is O-V-O with both the oxygen atoms in the leg 
direction. All angle values are listed in degrees. 



characteristic 


CaV 2 5 


MgV 2 5 


a, b, c 


11.35, 3.60, 4.89 


11.02, 3.69, 9.97 


dnn, d run g, d# e g(=6), ddiag 


3.03, 3.49, 3.60, 5.02 


2.98, 3.57, 3.69, 5.00 


a, /3, 7 


132.91, 102.94, 135.29 


117.57, 109.23, 141.15 
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where | inlma) (i denotes the site, n the main quan- 
tum number, I- orbital quantum number, m- magnetic 
number and a- spin index) are d-orbitals of transition 
metal ions. The density matrix is defined by: 
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where G a , 



, (E) = {inlma \ (E - H)' 1 \ inlm a) 
are the elements of the Green function matrix, N a = 
Tr(n° lm r) and N = + U and J are the screened 
Coulomb and exchange parameters. The U m m'm"m"' is 
the screened Coulomb interaction among the nl elec- 
trons which can be expressed via integrals over com- 
plex spherical harmonics and U, J parameters. For the 
Ca(Mg)V2 05 compounds the values of these parameters 
were calculated to be [7=3.6 eV and J=0.88 eV via the 
so-called "supercell" procedure [[l4| (in the "supercell" 
calculation only the xy-orbitals were considered to be lo- 
calized so that all other <i-orbitals could contribute to 
the screening). The calculation scheme was realized in 
the framework of the Linear Muffin-Tin Orbital (LMTO) 
method |J| based on the Stuttgart TBLMTO-47 com- 
puter code. 

Based on the Green function method, the inter-site ex- 
change couplings can be derived as the second derivative 
of the ground state energy with respect to the magnetic 
moment rotation angle Ilfi|l2]: 
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where the spin-dependent potentials I are expressed in 
terms of the potentials of Eq. (Q) , 
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The effective inter-sublattice susceptibilities are de- 
fined in terms of the LDA+U eigenfunctions ip as 
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The LDA+U method is the analogue of the Hartree- 
Fock (mean-field) approximation for a degenerate Hub- 
bard model ||. While in the multi-orbital case a mean- 
field approximation gives reasonably good estimates for 




FIG. 4. The angular distribution of the V d electron spin 
density for two V atoms belonging to the rung. Oxygen pyra- 
mids enclosing the V atoms are also shown in the Figure. 

the total energy, for the non-degenerate Hubbard model 
it is known to underestimate the triplet-singlet energy 
difference (and thus the value of the effective exchange 
coupling Jij) by a factor of two for a two-site problem 

(Ehf = and E exact = where t <C U is inter-site 
hopping parameter). 

As we discussed in section II the main building block 
of the crystal structures of the Ca(Mg)V205 compounds 
are the V ions roughly in the center of a pyramid of oxy- 
gen ions. The relevant point group symmetry is C^ v . The 
five <i-orbitals of the vanadium ion transform according to 
the following irreducible representations: 3z 2 — r 2 (Ai), 
x 2 — y 2 (Bi), xy (B2) and (xz,yz) (E). The lowest en- 
ergy orbital is the V3<i-orbital of xy-symmetry (using a 
convention where the axes of the coordinate system is 
directed towards the oxygen ions), which is the orbital 
whose lobes point at the directions, where the overlap 
with the oxygen is the smallest as can be seen in Fig. 4. 

Due to the crystal field splitting, the degeneracy of the 
V3<i-shell is lifted and the single d electron of V 4+ ion oc- 
cupies this rry-orbital, which reminds us of the cuprates, 
with a single hole in the x 2 — y 2 -orbital. The important 
difference is that while in cuprates all copper atoms are in 
the same (x, y)-plane as the x 2 — y 2 -orbital, in these vana- 
dates the vertices of the pyramids point up and down al- 
ternately with respect to the basal plane. Thus the V ions 
in their centers are correspondingly above and below the 
central plane, as can be seen in Fig. 2. As the xy-orbitals 
are parallel to this plane, so the overlap (and hence the 
exchange couplings) are expected to be stronger for vana- 
dium ions situated on the same side of the ladder plane. 
We will show that this is indeed the case. In addition 
to this alternation, a tilting of the pyramids is present in 
the crystal structure of these compounds, which we shall 
see seriously influences the interactions. 

Another important difference with the cuprates is that 

TABLE II. Calculated exchange coupling parameters 
(in K). The "Minus" sign indicates ferromagnetic exchange. 





CaV 2 5 


MgV 2 5 


Jl 


-28 


60 


J2 


608 


92 


J3 


122 


144 


J4 


20 


19 
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the zy-orbital has a 7r-overlap with the in-plane oxygen 
atoms in contrast to a much stronger c-overlap in the case 
of Cu 2+ . Consequently, one can expect much weaker ex- 
change interaction in vanadates as compared to cuprates. 
However it is surprising that the spin gap in CaA^Os 
(616 K H), is larger than the typical values for the sim- 
ilar cuprate ladders («460 K Q). 

So in the problem under consideration there are two 
types of contributions to the exchange interaction pa- 
rameters Jij . The first one is due to the xy — xy orbital 
hopping, and as only this orbital is half-filled this contri- 
bution directly corresponds to the non-degenerate Hub- 
bard model and its value must be multiplied by a factor 
of two to correct the Hartree-Fock value. Other contribu- 
tions are due to the hoppings to all other orbitals and as 
the mean-field approximation is much better for multi- 
orbital model this part can be used without modification. 

As mentioned earlier, the strongest interaction must 
be between V atoms which are situated on the same side 
of the plane (above or below) [see Fig. 1]. These atoms 
form ladders with interactions along the rung and the leg 
of the ladder denoted as J2 and J3 respectively and the 
interaction between the ladders as J\ (the notations are 
chosen to reflect the inter-atomic distances; the shortest 
one is between the atoms on different sides of the plane). 

Our calculated values of the exchange couplings are 
presented in Table II. It can be immediately seen that 
indeed the strongest interactions are between atoms on 
the same side of the plane (the ladder exchanges J2, 
J3). There is very strong anisotropy between the ex- 
change interactions along the rung ( J2=608 K) and the 
leg (J3=122 K) for CaV 2 5 . However for MgV 2 5 the 
rung ( J2=92 K) and the leg ( J3=144 K) exchange inter- 
action parameters are comparable in size. 

Our results suggests CaV205 is a system of weakly 
coupled dimers along the rung of the ladder with a very 
strong interaction inside the dimer. The analysis fl6j l 
based on fitting the results of model calculations to the 
experimental susceptibility measurements for CaV205 
confirms the coupled dimer picture and one of the ob- 
tained set of parameters (J2=665 K, J3=135 K, Jl=- 
25 K) is very close to our ab-initio calculated parameters 
values. However for MgV205 our calculations suggest 
jf=1.53 and =2.40 which puts MgV 2 5 outside the 
scope of the ladder limit, consistent with the helical or- 
dered gapless phase according to the phase diagram ob- 
tained by the Schwinger-boson mean field theory [p"7| . 
Recently jlq] the exchange parameters for CaV205 and 
MgV2 05 obtained in the LDA+U method were used 
for the calculations of the uniform susceptibility of the 
Heisenberg model by the quantum Monte Carlo method. 
The results agree very well with the experimental mea- 
surements, and particularly very good agreement has 
been found for CaX^Os. 



IV. LDA BAND STRUCTURES AND HOPPING 
INTEGRALS. 



In this section, we shall investigate the origin of the 
strong anisotropy of the rung and leg exchange interac- 
tions in CaV 2 05 and its absence in MgV205 using LDA 
band structure calculations. We shall use a systematic 
downfolding scheme to obtain an effective single (or few) 
band model Hamiltonian capable of reproducing the de- 
tails of the LDA bands close to the Fermi level. We shall 
extract the various hopping integrals which in turn could 
be related to the exchange interactions that we have cal- 
culated in the preceding section. 

Fig. 5(a) and Fig. 5(b) shows the energy bands and 
density of states (DOS) respectively for CaA^Os. The 
bands are plotted along the various high symmetry 
points jl9| of the Brillouin zone corresponding to the 
primitive orthorhombic lattice. Similarly in Fig. 6(a) and 
Fig. 6(b) we show the energy bands and density of states 
respectively for MgV2 05. The bands are now plotted 
along the various directions of the Brillouin zone cor- 
responding to base centered orthorhombic lattice. All 
the energies in the figures are measured with respect to 
the Fermi level of the respective compounds. In both 
compounds the bands below -3 eV have predominantly 
oxygen 2p character and are separated from the V d com- 
plex by a gap. From -1 cV to 3 eV the bands with V 3d 
characters are spread. In Fig. 5(c) and Fig. 6(c) we show 
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FIG. 5. Top left (a) TBLMTO-ASA energy bands; Top 
right (b) Density of states; bottom left (c) xy orbital pro- 
jected band structure and bottom right (d) Band structure of 
the effective four band model for CaV^Os. 
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FIG. 6. Top left (a) TBLMTO-ASA energy bands; Top 
right (b) Density of states; bottom left (c) xy orbital pro- 
jected band structure and bottom right (d) Band structure of 
the effective four band model for MgV^Os. 

the band structure of CaV2 0s and MgV2 0s respectively 
but projecting out the V xy orbital character, so that 
the fatness in each figure is proportional to the character 
of the orbital V xy (where x runs along the rung and y 
runs along the leg of the two dimensional ladder system) 
in the wavefunction. The orbital analysis of the bands 
or the so called fat bands confirm that the four lowest 
bands of the V 3d manifold are predominantly formed 
by V xy orbitals, consistent with the crystal field argu- 
ments presented in the previous section. The crystal field 
splitting between the xy-orbital and other 3<i-orbitals is 
so strong, that in the LDA band structure, particularly 
for MgV205, the xy-bands are separated from the rest 
of the V 3d-bands by a small energy gap as can be seen 
in Fig. 5b. With one subset of energy bands so well sep- 
arated from the rest, one can hope that a tight binding 
model with a single xy-orbital per V site (i.e. four band 
model because there are four V atoms in the unit cell) 
should provide a good approximation to the full band 
structure close to the Fermi level. In order to achieve this, 
the third generation TB-LMTO downfolding method |t]] 
has been employed to obtain a four band effective V-V 
model. The advantage of the new LMTO method over 
the previous ones is that the resulting downfolded Hamil- 
tonian is obtained in an orthogonal basis and accurately 
reproduces the full LDA bands close to the Fermi level. 
The energy bands obtained from the effective four band 
V-V model are shown in Fig. 5(d) and Fig. 6(d) respec- 
tively for CaV2C>5 and MgVzOs. It can be seen that the 



agreement with the LDA bands in Fig. 5(a) and Fig. 6(a) 
is remarkable. At this point it may be remarked that the 
downfolded orbitals are not thrown away but are included 
in the tails of the active LMTO's which are retained in 
the basis, i.e. in the tails of the V xy orbitals. In contrast 
to the fitting procedure often used to obtain tight-binding 
Hamiltonians the present method is deterministic and is 
free from adjustable parameters and also provides the 
information about the wavefunctions. The Fourier trans- 
form of the downfolded Hamiltonian H(k)^H(R) gives 
the effective hopping parameters for both compounds. 
Such an effective Hamiltonian is long ranged and has 
been called the physical Hamiltonian. We list in Table 3, 
all the hopping integrals which are relevant to understand 
the various exchange couplings presented in section III. 
It can be seen from Table 3, that analogously to the ex- 
change couplings the rung and leg hoppings of CaV20s 
are highly anisotropic while all the hoppings are of com- 
parable strength in MgV205. 

The exchange interaction parameter for the Hubbard 
model with strong Coulomb interaction can be estimated 
as J = 4t 2 /U, where t is the hopping parameter and U 
is the Coulomb interaction. For CaV20s the hopping 
along the rung of the ladder is 0.252 eV and £7=3.6 eV, 
which gives J=816 K (our LDA+U calculation gives 
J=608 K). The ratio of the rung and leg exchange pa- 
rameters (J2/J3) is equal to 4.98 while the ratio of the 
square of hopping parameters calculated in the down- 
folding procedure is equal to 6.25, confirming the strong 
anisotropy of the exchange couplings in CaV2C>5. The 
rung and leg hopping parameters for MgV205 are nearly 
equal, which again agrees with the LDA+U estimate for 
the exchange couplings in MgV20s. 

The hopping integrals extracted from the effective four 
band V-V model are consistent with the exchange cou- 
plings calculated from the LDA+U method in the pre- 
vious section. We shall now employ the downfolding 
method to explore the reason for the so different exchange 
as well as hopping integrals along the rung and leg for 
CaV205 in comparison to MgV20s although for both 
compounds the vanadium oxide planes have nearly the 
same geometry. 

It has been argued [^0| that at least the rung J2 and 
the leg J3 exchange integrals are mediated by the super- 
exchange mechanism through the O 2p orbitals and the 
size of this exchange integral primarily depends on the 



TABLE III. Calculated hopping parameters in an effective 
four band V-V model (in eV). The notations for the hopping 
parameters are same as that for the exchange interactions. 





CaV 2 5 


MgV 2 5 


tl 


0.076 


0.128 


t2 


0.252 


0.114 


t3 


0.101 


0.109 


a 


0.056 


0.069 
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hopping integral t p a between the xy orbital of the vana- 
dium and the p orbitals of the oxygen. According to the 
canonical band theory the structural difference between 
the two compounds should account for the observed dif- 
ferences. It should be noted that such super-exchange 
processes are explicitly taken into account in our effec- 
tive V-V model in the process of downfolding. 

Moreover in contrast to the cuprates, where the effec- 
tive hopping predominantly originates from the c-overlap 
of the Cu3d-orbitals with the oxygen 2p orbitals, for the 
vanadates the it overlap with the oxygen orbitals and the 
direct V 3d — 3d overlap could be of the same order of 
magnitude. 

The preceeding discussion suggests a model with oxy- 
gen pa; orbital along the leg and oxygen p y orbital along 
the rung in addition to the V xy orbitals should be a 
good starting point to understand these materials. Such 
a model Hamiltonian with more orbitals is usually short 
ranged and will be referred to as a Chemical Hamilto- 
nian, as it is expected to possess the necessary degrees of 
freedom so that its tight binding parameters behave in a 
meaningful way when the structure is deformed and when 
we proceed to study similar materials. Accordingly, we 
have extracted all the hoppings in a tight-binding Hamil- 
tonian where in the basis we have retained only oxygen p x 
orbitals along the leg, oxygen p y orbitals along the rung 
and all the V xy orbitals. All other orbitals were down- 
folded. Our calculation shows Y xy -O p hoppings (t p( j) 
are consistent with the prediction of the canonical band 
theory. In fact the ratio t r p ™ 19 for CaVzOs to MgV205 
is given to be 1.21 according to canonical band theory 
while our downfolding method yields 1.11. Further in this 
model the anisotropy between the rung and leg hopping is 
absent in CaV205. However for MgV20s the direct V xy - 
Y xy hopping along the rung is found to be very small 
in comparison to CaV205. This is not consistent with 
canonical band theory as the V-V distances in CaV2 05 
are smaller in comparison to MgVzOs (see table 1). In 
order to overcome this problem we tried to include more 
orbitals in the basis, as has been done for the High T c 
cuprates ID and the Ladder Cuprate SrCu 2 3 ||, how- 
ever the direct V-V as well as V-0 hoppings remained 
nearly the same as in the V-0 model discussed above. 
The reason a chemical Hamiltonian could not be defined 
for the vanadates in the same footing as the cuprates may 
be attributed to the complicated geometry of the vana- 
dates. As a consequence the orbitals are deformed in the 
process of downfolding thereby ruling out the validity of 
the simple canonical band theory 

However in this paper we have adopted the following 
strategy to overcome this problem. The chemical in- 
tuition suggests that the anisotropy of the leg and the 



rung exchange interactions in CaV20s and its absence 
in MgV2 05 may be attributed to the following: (a) the 
chemical composition of the compounds particularly the 
smaller ionic radii of Mg in comparison to Ca, (b) differ- 
ence in the crystal structure of these materials. 

In order to explore these effects we have considered 
three different models for CaV2 05 and have calculated 
the exchange interactions by the LDA+U method as ex- 
plained earlier and all the hopping parameters for the 
effective four band V xy -V xy model. The first model, re- 
ferred to as model 1, Ca is replaced with Mg in CaV2 0s, 
to examine whether the chemical composition plays any 
role in determining the anisotropy of the exchange inter- 
actions as well as hoppings in these materials. In order to 
explore the role of crystal structure we have considered 
the following two models. Model 2 is same as model 1 
except now the V-V and V-0 distances are changed so 
that they are equivalent to MgV2 0s. Finally in model 3 
we have not only changed the V-V and V-0 distances 
but also the V-O-V angles are changed so that they are 
the same as in MgV205. The results of our calculation 
for the exchange couplings as well as hoppings for the 
effective V-V model are summarized in Table 4. 

From Table 4 we conclude the following: The calcu- 
lation of the exchange coupling and the hopping param- 
eters in model 1 suggests that the change in chemical 
composition, i.e replacing Ca with Mg do not influence 
the leg and rung anisotropy as seen in CaV2 0s. Sim- 
ilarly from model 2 we conclude that bond lengths do 
not play any role in deciding the observed anisotropy 
between the leg and rung exchange interactions as well 
as hoppings in CaV20s. However calculations based on 
model 3 clearly shows as soon as the V-O-V angles are 
changed the exchange couplings as well as the effective 
hopping parameters are influenced appreciably. In this 
case, the rung exchange coupling and also the bare hop- 
ping is even smaller in comparison to the leg. We obvi- 
ously recover the values obtained for MgV20s as soon as 
the primitive orthorhombic stacking is changed to base 
centered orthorhombic stacking. These calculations sug- 
gests that the different tilting angle of the VO5 pyramids 
is the cause for the strikingly different magnetic behavior 
of the two vanadates considered here. 



TABLE IV. Calculated exchange couplings [in K] and hop- 
ping parameters [in eV] for the models as described in the 
text. 



System 


J2 


J3 


t2 


t3 


Modell 


320 


92 


0.169 


0.065 


Model2 


466 


57 


0.199 


0.045 


Model3 


24 


143 


0.053 


0.107 
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V. CONCLUSIONS 

We have used the LDA+U method to compute the 
exchange couplings in the layered vanadate compounds 
CaV2C>5 and MgV20s. Our calculation shows that a 
strong anisotropy exists between the rung and leg ex- 
change couplings for CaV20s thus making it a system 
of weakly coupled dimers along the rung with strong in- 
teraction inside the dimer, characterized by a large spin 
gap. On the other hand the rung and leg exchange 
couplings are found to be of comparable strength for 
MgV205 making it a small spin gap system. We have 
applied the recently developed third-generation LMTO 
downfolding method and the subsequently Fourier trans- 
formed the downfolded Hamiltonian to extract the tight- 
binding hopping parameters between effective V xy -V xy 
orbitals for CaV205 and MgV205, as well as for three 
different model systems. We conclude, that the stronger 
tilting of the VO5 pyramids in MgV205 crystal structure 
in comparison to CaV20s is the reason that the exchange 
interactions along rung and leg are nearly identical in 
MgV205 while they are anisotropic in CaV205 leading 
to the strikingly different magnetic properties of these 
materials. 
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